Passing current through touching molecules 
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The charge flow from a single Ceo molecule to another one has been probed. The conformation 
and electronic states of both molecules on the contacting electrodes have been characterized using 
a cryogenic scanning tunneling microscope. While the contact conductance of a single molecule 
between two Cu electrodes can vary up to a factor of three depending on electrode geometry, the 
conductance of the Ceo-Ceo contact is consistently lower by two orders of magnitude. First-principles 
transport calculations reproduce the experimental results, allow a determination of the actual Ceo- 
Ceo distances, and identify the essential role of the intermolecular link in bi- and trimolecular chains. 

PACS numbers: 73.63.-b, 68.37.Ef, 61.48.-c 



Intermolecular charge transport is central to numerous 
research fields. In biology electron hopping and tunneling 
processes between molecules play a vital role.^'^ More- 
over, tunneling processes between molecular materials 
have opened new perspectives towards the realization of 
efficient molecular sensors and solar cells In a parallel 
direction the conductance properties of point contacts)^ 
single atoms^ or single molecules^ are intensely being in- 
vestigated, and give a detailed view of charge transport 
through individual nanoscopic objects. Recently, experi- 
ments realized on ID extended moleculeSj^ single conju- 
gated polymers^ and DNA wiresi^ have been reported. A 
critical issue is now to understand and control the charge 
transfer from a single molecule to another one. 

Here we probe the electron current passing a chain of 
two Cgo molecules suspended in a STM junction, where 
the orientation and electronic states of both molecules 
have been characterized before connecting them with 
atomic-scale precision. The experimental results are 
complemented by first-principles transport calculations 
which give access to the distance-dependent nature of 
the intermolecular electron transport and predict the evo- 
lution of the transport properties with molecular chain 
length. 

The experiments were performed with a low- 
temperature STM operated at 5.2 K in ultrahigh vac- 
uum (below 10~* Pa). Au(lll) and Cu(lll) samples 
and etched W tips were prepared by Ar+ bombardment 
and annealing. As a final preparation, W tips were in- 
dented into the sample surface to coat them with surface 
material. Cgo molecules were deposited from a Ta cru- 
cible at a rate of ~ 1 ML/min as monitored by a quartz 
microbalance. During deposition a residual gas pressure 
in the 10^^ Pa range was maintained and the sample was 
kept at room temperature. The data shown correspond 
to a coverage of approximately 0.2 Cgo monolayers. All 
images were recorded in a constant-current mode. 

Increased image resolution with molecule-covered STM 
tips has repeatedly been reported . ^^d^d^ jjowever, no 
detailed information about the molecular orientation or 
their electronic properties was available. To realize con- 



FIG. 1: STM images (/ = 10 nA; V = 2.5 V; 14.3 x 10.5 nm^) 
of a Au(lll) surface partially covered with Ceo molecules 
(lower right) obtained with a (a) metal and (b) Ceo tip over 
the same area. Gold adatoms (a) and a small gold cluster (/3) 
of two or three adatoms are discernible. 



trolled molecular contacts, these details are decisive. 
We therefore used Ceo molecules as their orientation 
can be determined from sub-molecularly resolved STM 
iui'AgesM^^ Figure[T^ shows a STM image, recorded with 
a metallic tip, of an array of Cgo on a Au(lll) surface. 
Two Ceo orientations are observed, which are typical of 
a (2^/3 x 2V3)R30° Ceo superlattice^iS 

To attach a Ceo molecule to the tip, the metallic tip 
was placed over a target molecule and the sample voltage 
was varied from 2 V to 0.01 V and back at a constant 
current / — 100 nA. The success of this procedure can be 
verified from the removal of the molecule from the sub- 
strate (e.g., missing Ceo in Fig. [TJd indicated by a black 
arrow) . To further characterize Ceo tips, structures com- 
posed of one (a) and two or three {(3) Au adatoms had 
been deposited by slight contacts of the metallic STM tip 
with a clean surface area (Fig. [1^, upper left corner) 
The image of Fig. [T]3 was obtained with a Ceo function- 
alized tip over the same area. The Au clusters, which 
appear round and featureless with a metallic tip, exhibit 
a complex pattern which matches the highest occupied 
molecular orbital (HOMO) of Ceo-^"^ Obviously the Au 
adatoms work as tips for "reverse" imaging of Ceo 
the tip and provide direct access to the orientation of the 
molecule, e. g., in Fig.[T]D, a 6:6 bond of the Ceo tip is fac- 
ing the surface. While this technique has previously been 
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FIG. 2: DifFerential conductance {dl/dV) spectra acquired 
over (a) a Ceo molecule with a metallic tip and (b) the bare 
metal with a Ceo tip (Set point I — 10 nA; V — 2.5 V; rms 
modulation 10 mV). (c), (d) and (e) are reverse STM images 
(/ = 10 uA; 5 X 4.2 nm'^) acquired over a and /3 type atomic 
clusters with a Ceo tip at the voltages corresponding to peaks 
in the spectra of (b). (f), (g), and (h) are conductance maps 
acquired simultaneously with each picture (rms modulation 
50 mV). Reverse STM images (i) to (1) of a Au adatoms with 
a Ceo tip (/ = 10 nA; V = 2.5 V; 1.5 x 1.4 nm^). Between 
image acquisitions, the was approached to the surface so as 
to reach / ~ 1/iA, where reorientation of the molecule at the 
tip occurred. 



used to determine the number of molecules adsorbed on 
a STM tip,^^ the characterization of submolecular struc- 
tures was not reported. 

To monitor the density of states of Cgo tips, Fig. [5] dis- 
plays conductance spectra obtained with (a) a metallic 
tip on Ceo and (b) a " reverse" spectrum recorded with a 
Ceo tip on bare Au. The spectral peaks are characteris- 
tics of the molecular orbitals of Ceo on Au(lll)ji^ The 
spectra are almost perfect mirror images of each other re- 
flecting that the electronic state of Ceo at the tip is closely 
related to those of Ceo on the surface. "Reverse" images 
and conductance maps of atomic sized clusters exhibit 
submolecular patterns (Fig. [2t to h) which are typical 
of the lowest unoccupied molecular orbitals (LUMO-f 1, 
V w -2.5 V and LUMO, F « -1 V) and HOMO {V ^ 2 
y-j^jie Once a molecule is attached to the tip it is pos- 
sible to change its orientation by passing current of up 
to « 1/iA as demonstrated in Figs. [2j to 1. The molec- 
ular patterns obtained correspond to different Ceo ori- 
entations at the tip. While this sequence demonstrates 
control over the orientation of the tip molecule it also 
highlights an instability of these tips at high currents, 
which, therefore, were not suitable for the intended con- 
tact experiments. We repeated the previous experiments 
on Ceo deposited on a Cu(lll) substrate (Fig. [3^) where 
the binding of Ceo is stronger Transfer of Ceo from the 
Cu(lll) to the Cu-covered tip remains feasible, although 
the procedure is less reproducible than for Au(lll) and 
/iA-currents are required. As in the Au case, the struc- 



tural and electronic properties of the Ceo-tips have been 
characterized (Fig.[3]D). 

Following their characterization, metal and Ceo 
tips were approached to Ceo molecules and pristine 
Cu(lll) areas and conductance-distance [G(2;)] data 
were recorded (Fig.[3J;). No reorientation of the molecules 
occurred. Curve 1 was obtained with a sharp metallic tip 
approaching a Ceo molecule. The right part of the trace 
corresponds to the tunneling range. Contact is indicated 
by an inflection of the trace, which deflnes a contact con- 
ductance of 0.3 Go (Go — 2e'^ /h is the conductance quan- 
tum) in agreement with previous measurements on sim- 
ilar systems Curve 2 represents a measurement with 
a Ceo tip approaching pristine Cu. Surprisingly, the con- 
tact conductance of 1.0 Go is substantially higher than 
with Ceo on the surface. A conductance of « 1.5 Go 
at the contact constituted an upper limit which was not 
exceeded for different Ceo tip orientations. 

To understand the measured conductance traces flrst- 
principles transport simulations were carried out for rep- 
resentative junctions containing Ceo molecules between 
Cu(lll) electrodes. We modelled the fuUerene junctions 
by supercells with one or more Ceo molecules bridging 
a 4 X 4 representation of a slab containing 13 Cu(lll) 
layers. The electronic structure was determined with 
the SIESTA pseudopotential density functional theory 
(DFT) code^^ to calculate the transport properties for 
the TranSIESTA setup.^^ (For details cf. Ref. Case 
1 was modeled with a Ceo adsorbed on a hexagon on 
the substrate side centered underneath a Cu adatom on 
the tip side (cf. inset to Fig. [St), and case 2 with a 
Ceo adsorbed on a 5:6 bond on the tip side facing a 
clean Cu(lll) surface. For a transparent interpretation 
of the experiments we have considered the tip-molecule 
(sample-molecule) separation as the only variable in case 
1 (case 2), and full geometry relaxations were not per- 
formed. Except for the structural rearrangements ex- 
pected with a sharp metallic tip (case l)r^ this approach 
reproduces and explains the observed traces. The calcu- 
lated zero-bias conductances (Fig. [3j:) enable a calibra- 
tion of the absolute distances z (outermost Cu atom to 
Ceo-center along the surface normal) between tip (sam- 
ple) and molecule in case 1(2) by aligning the tunnel- 
ing part of the traces. Comparison of cases 1 and 2 
shows that for a given distance z, depending on the ge- 
ometry of the molecule-electrode interface, the conduc- 
tance of a single Ceo junction can vary by a factor of 
3 (10) under contact (tunneling) conditions. The con- 
ductance of the Ceo/Cu(lll) junctions is dominated by 
the molecular LUMO resonances that lie closest to the 
Fermi energy Ep. The theoretical maximum is there- 
fore 3 Go corresponding to three fully open conductance 
eigenchannels.^^ Indeed, a decomposition T — Ti of 
the total transmission T = T{Ep) into eigenchannel con- 
tributions {Ti\ confirms that the three most transmit- 
ting channels carry about an order of magnitude more 
current than the fourth. For the sharp-tip contact (case 
1 in Fig. 3) the transmissions in contact are of the order 
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FIG. 3: (a) dl/dV spectra acquired over (a) a Ceo molecule 
on Cu(lll) with a metallic tip and (b) the bare metal with 
a Ceo tip (Set point / = 10 nA; 1/ = 2 V; rms modulation 
10 mV. The insets show (a) a typical STM image (/ = 10 
nA; 1/ = 2 V; 2 X 1.6 nm^) of a Ceo array on Cu(lll), where 
all molecules expose hexagons to vacuum, and (b) a reverse 
STM image (7 = 10 nA; V = -2Y; 2.7 x 1.8 nm^) of the Ceo 
tip used for the contact experiment (a 5:6 bond is exposed to 
the surface). The inset to (c) displays sketches of the contact 
experiments performed by approaching (1) a sharp metallic 
tip to a Ceo adsorbed on a hexagon on Cu(lll), a 5:6 oriented 
Ceo tip (2) to the bare Cu(lll) surface and (3) to a Ceo 
adsorbed on a hexagon, (c) Experimental (solid lines) and 
calculated (symbols) conductances of the junctions versus the 
tip-molecule distance (1) (Set point / = 0.1 nA; V = 100 mV) 
and (2) (Set point 7 = 10 nA; = -100 mV) and the Ceo 
to Ceo center distance (3) (Set point 7 — 0.1 nA; V = 200 
mV). The calculated repulsive force (crosses) between two Ceo 
molecules suggests an elastic deformation of the junction at 
small separations that maps real molecule-molecule distances 
(open triangles) with apparent distances (filled triangles), see 
text. 



{T,} « {0.12,0.08,0.04,0.004}, hence the majority of an 
incoming electron wave is being reflected in this type of 
junction. Contrary, for the Cgo-tip contact (case 2 in Fig. 
3), three channels are much more open, theoretically in 
one case as much as {TJ « {0.97, 0.87, 0.57, 0.02}. 

To find out where the electrons are being scattered. 
Fig. visualizes the most transmitting eigenchannel 
wave function for the dilferent contacts.-- The isosur- 
faces are colored according to the phase as explained 
in the figure caption. Since the absolute square of the 
wave function corresponds to the density of the travers- 
ing electrons, the magnitude of the lobes gives an idea 
where the electron wave travels. In case 1 (sharp tip) 
the current is scattered at the single-atom contact to the 
molecule as indicated by the standing wave pattern at 
the tip side. In case 2 (Ceo-tip) the channel is almost 
perfectly open and the wave is propagating with essen- 




FIG. 4: (a) Visualizations of the most transmitting eigenchan- 
nel wave function (incoming from above) around the Fermi 
energy Ef- The isosurfaces of the real (imaginary) part of 
the wave function (with sign) are colored in white/dark blue 
(orange/light blue), (b) Calculated transmission functions for 
suspended molecular chains made of one, two, or three Ceo 
molecules. The molecular orientations correspond to (1) an 
adatom vs. a hexagon, (2) a 5:6 bond vs. a flat surface, (3) 
a hexagon vs. a 5:6 bond, and (4) a hexagon vs. a 5:6 bond 
vs. a hexagon. 



tially equal amplitude on either side of the molecule. To 
disentangle the effects of different molecular orientation 
as well as of different atomic contacts on the conduc- 
tance, we have carried out separate calculations with the 
hexagon orientation contacted with a flat tip. Specifi- 
cally, between flat Cu(lll) electrodes the 5:6 orientation 
was found to conduct slightly less than the hexagon ori- 
entation, e.g., at 2 = 6.68 A the conductance is 19 % 
lower for the 5:6 configuration. Moreover, experimental 
data from various Ceo orientations^ show that 0.3 Go 
(case 1) is already an upper limit of the conductance of 
Ceo on a Cu surface contacted with a sharp STM tip. We 
therefore conclude that the higher conductance of 1.0 Go 
(case 2) is due to the multiple-atomic contact which en- 
sures a better connection between the molecule and the 
electrode. This characterization of the metal-molecule 
contact could be valuable for fuUerene-based anchoring 
strategies for molecular electronics. 

Finally, the Ceo tip of Curve 2 in Fig. ^ was also ap- 
proached to a Ceo molecule on the substrate. Yazdani et 
al. used a similar method to measure the conductance 
of diatomic Xenon chain. Employing the molecular ori- 
entations determined in the experiment, this case was 
modeled with a Ceo adsorbed on a hexagon on the sub- 
strate side under a Ceo adsorbed on a 5:6 bond on the 
tip side. Note that the displacement axis now shows the 
Ceo to Cgo center distance. The observed conductance 
trace (Curve 3 in Fig. [5J;) varies smoothly from tunnel- 
ing to contact at a molecular separation of « 1 nm. The 
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contact conductance of « 0.01 Go is an order of magni- 
tude smaller than expected for a Cgo dimeri^l In contrast 
to the experimental observation of a plateau, our model 
predicts an exponential dependence of the conductance 
on the molecule-molecule separation (open triangles in 
Fig. and no significant infiuence of the molecule- 
surface distance (see Ref . ITol) . This difference is due to 
an intermolecular repulsion at small distances that de- 
forms the contact. Therefore, beyond the point of con- 
tact, the experimental data reflect apparent molecular 
separations, the actual distances being somewhat larger. 
To take into account the elastic deformation of the junc- 
tion, repulsive forces were estimated from our DFT cal- 
culations (crosses in Fig. [3t) . By renormalizing the theo- 
retical z-coordinates according to the compliance of two 
soft molecule-surface segments (effective elastic constant 
fccff = 7 eV/A),-i^ we obtain a good agreement with the 
experimental trace (filled triangles in Fig. (St) . While it is 
possible to further improve this agreement by consider- 
ing the elasticity of the other parts of the systemji^ the 
essential feature of the experiment is already captured 
by inclusion of just the softest segment. Interestingly, 
the onset of elastic deformation coincides with the inter- 
molecular distance of 1.004 nm in Cgo crystals, which is 
controlled by van der Waals bonding and electrostatic 
repulsion.— 

The picture emerging for chains of two Cgo molecules 
is that the transport processes are mainly sensitive to the 
molecule-molecule interface. It is further supported by 
Fig. [4^3 which shows isosurfaces of the dominant eigen- 
channel with little weight on the lower molecule. This 



is due to a reduced wave function amplitude beyond the 
Cgo^Cgo interface, which thus acts as conductance bot- 
tleneck. Within the chain the intermolecular distance is 
limited by electrostatic repulsion. In this way the experi- 
ment is probing how current passes through two touching 
molecules, the properties and the nature of both being 
controlled and tunable. 

Using the characteristic C6o~C6o contact distance de- 
termined above, the transport through a three-Ceo chain 
was calculated. In this case the dominant eigenchannel 
Fig. 2^4 is strongly attenuated along the chain as revealed 
by the absence of lobes on the lower molecule. The trans- 
mission functions of the molecular chains Fig. |4j3 reveal 
the opening of a ~ 1.5 eV gap around Ep, and hence pre- 
dict a rapid evolution towards the insulating character of 
an infinite Cgo chain. 

In summary, the contact conductance for single Cgo 
junctions can vary up to a factor of three depending on 
the molecule-metal interfaces, thus corroborating the no- 
tion of good and bad contacts. The current passing from 
one molecule to another one, however, is determined by 
the molecule-molecule interface. Our experimental ap- 
proach can be extended to a range of molecules to ad- 
dress the influence of the molecule-molecule interactions 
on intermolecular charge transport. Moreover, through 
detection of photons emitted in a STM junction, a suit- 
able fluorescent molecule attached to the STM tip might 
prove useful as optically active probej22 
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I. SUPPLEMENTARY MATERIAL 

Theoretical details. The electronic structure for 
the supercells shown in the inset of Fig. 3c was calcu- 
lated with the SIESTA^ pseudopotential density func- 
tional theory (DFT) method and the generalized gra- 
dient approximation^ for exchange-correlation (XC). A 
standard single-zeta plus polarization basis was employed 
for Cu (0.15 eV energy shift, Vc < 6.9 a.u.) and a long- 
ranged, double-zeta plus polarization basis for C (0.02 eV 
energy shift, < 5.5 a.u.). Real-space grid integrations 
were carried out using a 200 Ry energy cutoff. The 3D 
Brillouin zone was sampled with a 2 x 2 x 1 Monkhorst- 
Pack k-mesh. The lattice constant for the Cu crystal 
was set to 3.70 A. An initial structure consisting of a Ceo 
molecule adsorbed with a hexagon on an hep hollow-site 
on a slab containing 7 Cu(lll) layers, was fully relaxed 
until all forces on the molecule, the adatom placed on 
the reverse side of the slab, and the surface layer, were 
smaller than 0.02 eV/A. Based on this configuration and 
the molecular orientations derived from the experiment 
we constructed the supercells shown in Fig. 3c for several 
different electrode separations. These supercells were not 
relaxed; forces up to a few eV/A were thus tolerated. 

The electronic structure from SIESTA was used to 
calculate the transport properties for the TranSIESTA'^ 
setup. The zero-bias conductance G = GoT{Ef) was 
derived from the transmission function T{E) calculated 
by Green's function techniques involving the Kohn-Sham 
Hamiltonian in the scattering region and self-energies 
representing atomistic, semi-infinite electrodes obtained 
from separate calculations for bulk Cu(lll). The trans- 
mission was sampled over a 6 x 6 k||-mesh for the 2D 
Brillouin zone. To check the parameters described above, 
we carried out calculations for a representative structure 
with increased number of k-points for both the electronic 
structure part as well as for the transmission. The varia- 
tions in the zero-bias conductance were < 6%. The eigen- 
channel visualizations, shown in Fig. 4a, were calculated 
at the Fermi energy Ep for the F-point according to the 
scheme presented in Ref.|23l As seen in Fig. 4b the trans- 
mission functions vary smoothly over the energy scale of 
the experimentally applied voltages. This justifies the 
comparison with the calculated zero-bias conductances. 

Elastic response of the Ceo— Ceo junction.. The 

theoretical data in Fig. 3c show that the conductance 
of the Ceo-Ceo junction depends exponentially on the 
intermolecular separation, approximately as G{z) oc 

g~2.6A prove that the conductance does not de- 

pend significantly on the molecule-substrate distance, we 
varied the electrode separation by ±0.50 A for a fixed 
(realistic) Cgo-Ceo distance of z — 9.54 A, only to find 
the conductance changed by less than 9%. In compari- 
son, varying the intermolecular separation by the same 
amount, the conductance changes by a factor 3.7. The 
experimental observation of a plateau after contact for- 
mation is therefore related to an elastic deformation of 
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FIG. SI: Conductance traces and intermolecular forces for 
the Ceo-Ceo junction with different estimates of the elasticity 
of the system and the Ceo-Ceo repulsion (with and without 
van der Waals contributions). 

the junction, driven by a strong intermolecular repulsion 
between the two molecules. 

To estimate the Ceo-Ceo repulsion in our structures 
we determined the average force {fP - fP)/2, where 

i^i*-* is the z-component of the sum of atomic forces on 
molecule i from our DFT calculations. Taking the setup 
with the largest intermolecular separation as a point 
of reference, the evolution of intermolecular repulsion 
within DFT is represented in Figs. 3c andlSlI For an elas- 
tic system, characterized by the effective elastic constant 
kcs arising from a series of springs (l/fccfr = 1/^i): the 
STM piezo distances can thus be related to intermolec- 
ular distances via Zpjezo = -Winter — F/kcs- One example 
of such coordinate renormalization is shown in Fig. 3c 
where only the two soft molecule-surface segments were 
included (/coff = fcm-s/2). One can refine this descrip- 
tion by including Ceo compressions (k^aoi 108 eV/A^ 
for each Ceo) and the compliance of the metal surfaces 
(^surf ~ 30 eV/A^ for each electrode), which leads to 
kcs = 4.4 eV/A^ and an a very good agreement with the 
experimental data, cf. triangles in Fig. ISII 

While short-range interactions are generally well- 
described in DFT with standard local approximations 
for XC, the long-range van der Waals (vdW) interac- 
tions are not easily incorporated in DFT because of 
its nonlocality. Therefore, to get an estimate of the 
vdW contributions to the forces in the Ceo-Ceo con- 
tact (case 3 in Fig. 3), we applied an empirical correc- 
tion following Ref. 0. In this scheme the interaction en- 
ergy between each pair of atoms {i, j), separated by 
the distance R, is corrected by an additional attrac- 
tive energy the form C/^^^w ^ - f^{R)Cl' / ^ where 

fl{R) = 1 - exp[-0.00075i?V(ri + r^f] is a damp- 
ing function that removes contributions for small sep- 
arations. The parameter is a characteristic covalent 
radius of atom i. We used tq = 0.7245 A (correspond- 
ing to the average bond length in our description of Ceo) 
and rcu = 1.308 A (derived from the Cu lattice con- 
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stant). For the Cg coefficient we used the London formula 
Cq = ^aittjlilj / [li + Ij) with experimental values^ for 
the polarizability ai and the ionization potential U of 
atom i. The strategy outlined above enabled us to re- 
produce the 3D fee crystal of Ceo with a lattice constant 
of 10.02 A and a bulk modulus only about 12% lower 
than the experimental value. Applying the vdW correc- 
tion to the Cgo-Cgo contact the forces shown in Fig. [ST] 
are obtained. Since the intermolecular repulsion is only 



slightly reduced by the attractive vdW forces, the coordi- 
nate renormalization of the conductance trace is mostly 
sensitive to the elastic constant fceff and not to the force 
estimate. Furthermore, since it is generally argued that 
the London form shows a tendency to overestimate the 
vdW interaction [see Ref. ? and references herein], we 
conclude that vdW contributions play a minor role for 
understanding the experimental conductance trace. 
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